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ABSTRACT 

We study the prospects for extracting detailed statistical properties of the Sunyaev-Zel'dovich (SZ) 
effect associated with large scale structure using upcoming multifrequency CMB experiments. The great- 
est obstacle to detecting the large-angle signal is the confusion noise provided by the primary anisotropics 
themselves, and to a lesser degree galactic and extragalactic foregrounds. We employ multifrequency 
subtraction techniques and the latest foregrounds models to determine the detection threshold for the 
Boomerang, MAP (several /iK) and Planck CMB (sub /iK) experiments. Calibrating a simplified biased- 
tracer model of the gas pressure off recent hydrodynamic simulations, we estimate the SZ power spectrum, 
skewness and bispectrum through analytic scalings and N-body simulations of the dark matter. We show 
that the Planck satellite should be able to measure the SZ effect with sufficient precision to determine 
its power spectrum and higher order correlations, e.g. the skewness and bispectrum. Planck should 
also be able to detect the cross correlation between the SZ and gravitational lensing effect in the CMB. 
Detection of these effects will help determine the properties of the as yet undetected gas, including the 
manner in which the gas pressure traces the dark matter. 

Subject headings: cosmic microwave background — cosmology: theory — large scale structure of 
universe 



1. INTRODUCTION 



It is by now well established that the precision measure 

ments of the cosmic microwave background expected from 
upcoming experiments, especially MAP and Planck satel- 
lite missions, will provide a gold mine of information about 
the early unive rse and the fun damental cosmological pa- 
rameters (e.g., Jungman ct al.| 1996). These experiments 
can in fact do so much more. With all-sky maps across the 
wide range of uncharted frequencies from 20GHz-900GHz, 
the secondary science from these missions will arguably be 
as interesting as the primary science. 

In this paper, we examine the prospects for extract- 
ing the large-scale properties of the hot intergalactic gas 
from multifrequency observations of the CMB. Inverse- 
Compton scattering of CMB photons by hot gas, kno wn as 
the Sunyaev-Zel'dovich (SZ; (Bunyaev k Zel'dovich| 1980) 
effect, leaves a characteristic distortion in the spectrum 
of the CMB, which fluctuates in the sky with the gas 
density and temperature. In the Rayleigh- Jeans (RJ) 
regime, it produces a constant decrement and with only 
low frequency measurements, the much larger primary 
anisotropics in the CMB itself obscure the flu ctuations on 
scales g reater than a few arcminutes (e.g., Goldberg & 
Sperge] 1999). The upscattering in frequency implies an 
increment at high frequencies and a null around 217GHz. 
This behavior provides a potential tool for the separation 
of SZ effect from other temperature anisotropy contribu- 
tors. 

Since both the SZ spectrum and the CMB spectrum are 
accurately known, one can expect that foreground removal 
techniques developed to isolate the primary anisotropies 
can be reversed to recover the SZ signal in the presence of 
noise from the primary anisotropies. Galactic and extra- 
galactic foregrounds will be more challenging to remove. 



He re we use the latest foreground models from Legmark ct 
al. (1999) that takes into account the fact that imperfect 
correlations in the foregrounds between frequency chan- 
nels inhibits our ability to remove them. Using foreground 
information together with the expected noise properties 
of individual experiments, one can determine the mini- 
mal detectable signal in each experiment and the upper 
limit achievable in the absence of detection. Experiments 
with sufficient signal-to-noise can extract precision mea- 
surements for the power spectrum and higher order statis- 
tics such as the skewness. Ultimately, they can provide 
detailed maps of the large-angle SZ effect. 

To assess the prospects for an actual detection, we must 
model the SZ signal itself. The SZ effect is now routinely 
image d in massive galaxy clusters (e.g., Carlstrom et &\ 
1996; Jones et a] 1993), where the temperature of the 
scattering medium can reach as high as 10 keV, producing 
temperature changes in the CMB of order 1 mK at RJ 
wavelengths. The possibility for detection of massive clus- 
ters in CMB satellit e data has already been discussed in 
sever al studies (e.g., |Aghanim ct ah , Hachnclt & Tegmark 



1996, Pointccouteau et a| 1998). Here, however, we are in- 



terested in the SZ effect produced by large-scale structure 
in the general intergalactic medium (IGM) where the gas 
is expected to be at <J IkeV in mild overdensities, leading 
to CMB contributions in the fjK range. 

It is now widely believed that at least ~ 50% of the 
present day baryons, when compared to the total baryon 
budget from big bang nucleosynthesis, are present in gas 
associated with hot large-scale structure which has re- 
mained undetec ted given its t empe rature and clust ering 
prop erties (e.g., Fukugita ct al 1998 ; pen fc Ostrikei 1999; 



Pen] 1999). Recently, jScharf ct al| (2000) has provided a 



tentative detection of X-ray emission from a large-scale 



2 



filament in one of the deep ROSAT PSPC fields; previ- 
ous a ttempts yielding u pper limits are described in Kulj 
& Bo frringerj (1999) and [Briel fc Henry) (1995). These re- 
sults are consistent with current predictions for the X-ray 
surface b rightne ss ba sed on numerical simulations (e.g., 
Ccn ct al 1995). Pen (1999) argued that non-gravitational 
heating of the gas to ~ IkeV is required to evade bounds 
from the soft X-ray background. These results suggest 
that the X-ray emission from this gas may be detectable in 
the near future with wide-field observations with Chandra 
X-ray Observatory^ and X-ray Multiple Mirror Mission^. 

On the theoretical front, hydrod ynamic simulat ions of 
the SZ effect co ntinue to improve ( da Silva et al. 1999; 
Refregier et~ai] 1999; [Beljak et al] 2000). As a consen- 



sus from these simulations of basic properties such as the 
opacity weighted gas temperature and average Compton 
distortion is still lacking, we will base our assessment of the 
detectability of the large-scale SZ effect on a simple param- 
eterization of the effect, based on a gas pressure bias model 
( Refregier ct al. 1999), crudely calibrated with the recent 
hydrodynamic simulations. We employ perturbation the- 
ory, non-linear scaling relations, and N-body simulations 
for the dark matter to assess the statistical properties of 
the signal. Properly calibrated, these techniques can com- 
plement hydrodynamic simulations by extending their dy- 
namic range and sampling volume. Currently, they should 
simply be taken as order of magnitude estimates of the 
potential signal. 

Throughout this paper, we will take an adiabatic cold 
dark matter (CDM) model as our fiducial cosmology. We 
assume cosmological parameters ft c — 0.30 for the cold 
dark matter density, fib = 0.05 for the baryon density, 
f^A = 0.65 for the cosmological constant, h = 0.65 for the 
dimensionless Hubble constant and a COBE- norm alized 
scale inva riant spectrum of primordial fluctuations ( Bunrj| 
& WHhc] 1997). 

The layout of the paper is as follows. In § H, we describe 
the foreground and primary anisotropy removal method 
and assess their efficacy for upcoming CMB experiments. 
In § ||, we detail the bias model for the SZ effect and calcu- 
late through perturbation theory, analytic approximations 
and numerical simulations, the low order statistics of the 
SZ effect: its power spectrum, skewness and bispectrum. 
In § I], having estimated the noise and the signal, we as- 
sess the prospects for measuring these low order statistics 
in upcoming experiments. We conclude in § |5| with a dis- 
cussion of our results. 

2. MODELING THE CMB AND FOREGROUND NOISE 

The main obstacle for the detection of the SZ effect from 
large-scale structure for angular scales above a few arcmin- 
utes is the CMB itself. Here the primary anisotropies dom- 
inate the SZ effect for frequencies near and below the peak 
in the CMB spectrum (see Fig. [I]). Fortunately, the known 
frequency dependence and statistical properties of primary 
anisotropies allows for e xtremely effect ive su btraction of 



their co ntribution (e.g., Hobson et al 1998; Bouchet & 



Gispert 1999). In particular, primary anisotropies obey 
Gaussian statistics and follow the blackbody spectrum pre 
cisely. 

1 http:/ /asc. harvard.edu 

2 http: / /astro. estec.esa.nl/XMM 



Table 1 



CMB Experimental Specifications 



1 / AIM . 1 




FWHM 


10 6 AT/T 


Boomerang 


90 


20 


7.4 




150 


12 


5.7 




240 


12 


10 




400 


12 


80 


MAP 


22 


56 


4.1 




30 


41 


5.7 




40 


28 


8.2 




60 


21 


11.0 




90 


13 


18.3 


Planck 


30 


33 


1.6 




44 


23 


2.1 




70 


14 


3.6 




100 


10 


4.3 




100 


10.7 


1.7 




113 


8.0 


2.0 




217 


5.5 


4.3 




353 


5.0 


14.4 




545 


5.0 


147 




857 


5.0 


6670 



NOTES. — Specifications used for Boomerang, MAP and Planck. 
Full width at half maxima (FWHM) of the beams are in arcminutcs 
and should be converted to radians for the noise formula. Boomerang 
covers a fraction ~ 2.6% of the sky, while we assume a usable frac- 
tion of 65% for MAP and Planck. In § W, in order to calculate the 
maximum signal-to-noise, we define a perfect experiment as one with 
no instrumental noise and full sky coverage. 



Perhaps more worrying are the galactic and extragalac- 
tic foregrounds, some of which are expected to to be at 
least comparable to the SZ signal in each frequency band. 
These foregrounds typically have spatial and/or temporal 
variations in their frequency dependence leading to imper- 
fect correlations between their contributions in different 
frequency bands. We attempt here to provide as realis- 
tic an estimate as possible of the prospects for CMB and 
foreground removal, given our incomplete understanding 
of the foregrounds. 

2.1. Foreground Model and Removal 



We use the "MID" foreground model of Tegmark et al 
(1999) and adapt the subtraction techniques found there 
for the purpose of extracting the SZ signal. The assumed 
level of the foreground signal in the power spectrum for 
three fiducial frequencies is shown in Fig. [j]. 

The foreground model is defined in terms of the covari- 
ance between the multipole moments at different frequency 
bands 3 



{a\t m ,{v')a f lm {v)) = d{{i/,v)8 w &« 



(1) 



in thermodynamic temperature units as set by the CMB 
blackbody. In this section, we will speak of the pri- 
mary anisotropies and detector noise simply as other fore- 
grounds with very special properties: 

3 A potential caveat for this type of modeling is that it assumes 
the foregrounds are statistically isotropic whereas we know that the 
presence of the Galaxy violates this assumption at least for the low 
order multipoles. We assume that 1 — f s ^ y ~ 0.35 of the sky is lost 
to this assumption even with an all-sky experiment. 
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Fig. 1. — Top: foreground contributions to temperature anisotropics (AT/T) 2 = 1(1 + l)Ci/2n from the various foregrounds (dust, free-free, 
synchrotron, radio and infrared point sources, and rotating dust) at three fiducial frequencies as labeled. The SZ signal (solid, unlabeled) is 
estimated with the simplified model of § 3. Bottom: residual foregrounds after multifrequency subtraction for Boomerang, MAP and Planck. 
The total includes detector noise and residual CMB. 



q noisc (^» = 8 In 20(f) V 



/AT 



(2) 



The FWHM= \/8 In 28 and noise specifications of the 
Boomerang, MAP and Planck frequency channels are 
given in Tab. 1. True foregrounds generally fall in between 
these extremes of perfect and no frequency correlation. 

The difference between extracting the SZ signal and the 
primary signal is simply that one performs the subtraction 
referenced to the SZ frequency dependence 



s(v) = 2 — — coth — , 
W 2 2 ' 



(3) 



where x = hv/kT cm b ~ ^/56.8GHz. Note that in the RJ 
limit s(y) — > 1 such that 



Cf>,i/) = *(i/)«(i/)Cf z 



(4) 



where Cf z is the SZ power spectrum in the RJ limit. 

Consider an arbitrary linear combination of the chan- 
nels, 

b = £ ~f—\ w { v i) a { v i) ■ ( 5 ) 

v V V 

Since the subtraction is done multipole by multipole, we 
have temporarily suppressed the multipole index. The co- 



variance of this quantity is 



g f (^j, Vj) 
s{Vi)s{Vj) 



(6) 

Minimizing the variance contributed by the foregrounds 
subject to the constraint that the SZ estimation be unbi- 
ased, we obtain 



const. 



(7) 



Defining the scaled foreground covariance matrix as 



(8) 



we solve for the weights that minimize the noise variance 

wocC^e, (9) 

where e is the vector of all ones e(ui) = 1. Finally we nor- 
malize the sum of the weights to unity = 1 to ob- 
tain an unbiased estimator. Our approach is same as min- 
imizing the foreground variance subject to the constraint 
that the recovered multipole is an unbiased estimate of 
the true SZ signal. As each channel is rescaled such that 
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Fig. 2. — Dependence of the residual noise rms on foreground 
assumptions expressed as a ratio to the fiducial model of Fig. [H. (a) 
Falsely assuming the foregrounds have perfect frequency coherence 
not only underpredicts the residual noise by a substantial factor but 
also leads to substantially more actual residual noise, (b) Multiply- 
ing the foreground amplitudes by 2 (power by 4) produces less than 
a factor of 2 increase in the residual noise. 



SZ signal corresponds to the RJ level, the weights sum to 
unity. 

The total residual noise variance in the map from the 
foregrounds per multipole is then 



Ni = w*C/w ; , 
and from each foreground component 
Nf = w*C[w, . 



(10) 



(11) 



Note that the residual noise in the map is independent 
of assumptions about the SZ signal including whether it 
is Gaussian or not. However if the foregrounds them- 
selves are non-Gaussian, then this technique only mini- 
mizes the variance and may not be optimal for recovery 
of non-Gaussian features in the SZ map itself. Bouchct 
et alQl995) have shown that similar techniques are quite 
effective even when confronted with non-Gaussian fore- 
grounds. This is a potential caveat especially for cases in 
which the residual noise is not dominated by the primary 
anisotropies or detector noise. We shall discuss methods 
to alleviate this concern in the next section. 

2.2. Detection Threshold 

The residual noise sets the detection threshold for the 
SZ effect for a given experiment. In Fig. [I], we show the 
rms of the residual noise after foreground subtraction for 
the Boomerang, MAP and Planck experiments as suming 



the "MID" foreground model from |Tegmark et al. (1999) 



With the Boomerang and Planck channels, elimination of 
the primary anisotropies is excellent up to the beam scale 
where detector noise dominates. As expected, the MAP 
channels, which are all on the RJ side of the spectrum, do 
not allow good elimination of the primary anisotropies. 

It is important not to assume that the foregrounds are 
perfectly correlated in fre quency, which is th e usual as- 
sumption in the literature (Hobson ct al 1998; Bouchet & 



Gispert 1999). There are two types of errors incurred by 
doing so. The first is that one underpredicts the amount 
of residual noise in the SZ map (see Fig. |J). The sec- 
ond is that if one calculates the optimal weights in equa- 
tion (^) based on this assumption the actual residual noise 
increases. For Planck it can actually increase the noise 
beyond the level in which it appears in the 100GHz maps 
with no foreground subtraction at all. The reason is that 
the cleaning algorithm then erroneously uses the contam- 
inated high and low frequency channels to subtract out 
the small foreground contamination in the central chan- 
nels. In Planck, the difference between the predicted and 
actual rms noise from falsely assuming perfect frequency 
coherence can be more than two orders of magnitude. 

For Boomerang and Planck, the largest residual noise 
component, aside from detector noise, is dust emission and 
is sufficiently large that one might worry that current un- 
certainties in our knowledge of the foreground model may 
affect the implications for the detection of the SZ effect. It 
is therefore important to explore variations on our fiducial 
foreground model. 

Multiplying the foreground rms amplitudes uniformly 
by a factor of 2 (and hence the power by a factor of 4), 
produces less than a factor of 2 increase in the residual 
noise rms as sh own in Fig. 0. Likewise, as discussed in 



Tegmark ct al. (1999), minor variations in the frequency 



coherence do not effect the residual noise much in spite 
of the fact that it is crucial not to assume perfect corre- 
lation. We conclude that uncertainties in the properties 
of currently known foregrounds are unlikely to change our 
conclusions qualitatively. There is however always the pos- 
sibility that some foreground that does not appear in the 
currently-measured frequency bands will affect our results. 

The fact that the residual dust contributions are com- 
parable to those of the detector noise for Boomerang and 
Planck is problematic for another reason. Since the al- 
gorithm minimizes to total residual variance, it attempts 
to keep these two main contributors roughly comparable. 
However the dust will clearly be non-Gaussian to some 
extent and one may prefer instead to trade more residual 
detector noise for dust contamination. One can modify the 
subtraction algorithm to account for this by artificially in- 
creasing the rms amplitude of the dust when calculating 
the weights in equation (^) while using the real ampli- 
tude in calculating the residual noise in equation ( |TT| ) . For 
example we have explored increasing the amplitude by a 
factor of 4 (power by 16) for the weights. The result is 
an almost negligible increase in total residual noise rms 
but an improvement in dust rejection by a factor of 3-4 
in rms. For Planck this brings the ratio of dust to total 
rms to ~ 10% and recall that the noise adds in quadrature 
so that the total dust contribution is really ~ 1% of the 
total. This more conservative approach is thus advisable 
but since it leaves the total residual noise rms essentially 
unchanged, we will adopt the minimum variance noise to 
estimate the detection threshold. 

Fig. H directly tells us the detection threshold per (I, m) 
multipole moment. Since the SZ signal is likely to have a 
smooth power spectrum in I, one can average over bands 
in I to beat down the residual noise. Assuming Gaussian- 
statistics, the residual noise variance 2Nf for the power 



5 



spectrum estimate is then given by 



band 



/sky 

'band 



i)jvr 



(12) 



where / s k y accounts for the reduction of the number of in- 
dependent modes due to the fraction of sky covered. The 
result for the three experiments is shown in Fig. ^. In the 
absence of a detection, they can be interpreted as the opti- 
mal 1 a upper limits on SZ bandpowers achievable by the 
experiment. Boomerang and MAP can place upper lim- 
its on the SZ signal in the interesting /iK regime whereas 
Planck can detect signals well below a /iK. 

This noise averaging procedure in principle implicitly as- 
sumes that the statistical properties of the residual noise, 
and by implication the full covariance matrix of the other 
foregrounds, is precisely known. In reality, they too must 
be estimated from the multifrequency data itself through 
either through the subtraction techniques discussed here 
or by direct m odeling; of the foregrounds in the maps. 
Tegmark et alj (1999) found that direct modeling of the 
foregrounds with hundreds of fitted parameters did not 
appreciably degrade our ability to extract the properties 
of the primary anisotropics. The main source of variance 
there was the cosmic variance of the primary anisotropics 
themselves whose properties are precisely known. Simi- 
larly here the main source of residual variance is either 
the primary anisotropies (for MAP) or detector noise (for 
Boomerang and Planck) and their statistical properties 
may safely be considered known. 

3. MODELING THE SZ SIGNAL 

In order to estimate how well the statistical properties 
of the SZ effect might be recovered with multifrequency 
CMB maps, we need to model the large-angle SZ effect 
itself. T he current stat e-of-th e-art in hydrody nami c simu- 



lation s (ia Silva et al. 1999; Refregier et al. 1999; Seljak 



et al. 2000) has reached a qualitative but not quantita- 



tive consensus on the statistical properties of the SZ cf- 
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Fig. 3. — Detection thresholds for the SZ effect. Error boxes 
represent the 1-tr rms residual noise in multipole bands and can be 
interpreted as the detection threshold. Also shown (dotted) is the 
level of the primary anisotropies that have been subtracted with the 
technique and the signal (dashed) expected in the simplified model 

of §y. 



feet. In addition, questions as to the heating of the gas 
from non-gravitati onal sources may even change the re- 
sults qualitatively (Pen 1999). Hydrodynamic simulations 
are also severely limited in the dynamic range and volume 
sampled. 

Given the current state of affairs, we believe it is use- 
ful to explore a parameterized model of the effect whose 
consequences are simple to calculate and which may be 
calibrated against hydrodynamic simulations as they con- 
tinue to improve. 

3.1. Bias Prescription 

In general, the SZ temperature fluctuation 9 = AT/T 
is given by the opacity weighted integrated pressure fluc- 
tuation along the line of sight: 



e sz (n,z/) = -2s{v) / dr 
Jo 



T7r(r, nr) . 



(13) 



r is the the comoving distance, r is the Thomson optical 
depth, overdots are derivatives with respect to r and the 
dimensionless electron pressure fluctuation is 



(14) 



One needs to model the statistical properties of ir, in par- 
ticular its power spectrum and bispectrum 



(7r(k)V(k')) 
(7r(k)7r(k')7r(k")> 



(27T) 3 (5 D (k- 
(27T) 3 (5 D (k- 



k')P.(fc), (15) 
■k' + k'O-B^fc, &',&"), 



as a function of lookback time or distance r. In principle 
we also need the unequal time correlators, but in practice 
these do not play a role as we shall see. 

By analogy to the familiar case of galaxy clustering, it 
is reasonable to suppose that the pressure fluctuations de- 
pend locally on the dark matter density and hence are bi- 
ased tracers of the dar k matter density in the linear regime 
( Goldberg fc Spergel 1999). Hence the statistical proper- 



ties follow from those of the dark matter distribution 



P„(k;r)nb v (r) 2 P 5 (k;r), 
B„(k, k', k"; r) a b„{r) 3 B s {k, k' , k"; r) 



(16) 



We have restored the time dependence since the bias will 
be time dependent even in the linear regime and must be 
extracted from simulations. In general, the bias parameter 
for the power spectrum and bispectrum need not be the 
same even in the linear regime since the bispectrum au- 



tomatic ally involves higher order corrections (Fry & Gaz- 
tanaga 1993). For estimation purposes here we will take 

them to be equal. 

Following |Goldberg fc Spergel (1999), we chose the form 



K{r) =&*(0)/(l + *), 



(17) 



as motivated by findings that the average gas temperature 
drops off roughly by this factor. We normalize the value 
of the bias parameter today by comparison with recent 
hydrodynamic simulations. It is conceptually useful to 
separate the bias into two factors: 



6.(0) = 



k B T e (0) 



m e c* 



bs . 



(18) 



6 



i.e. an opacity-weighted average temperature and a bias 
parameter for the gas density at that temperature. In 
Rcfregicr et al.| (1999), for our fiducial ACD M cosmology, 
the bias b$ was found to be ~ 8 to 9, while in Seljak et al. 
(2000) it was found to be in the range ~ of 3 to 4. In both 
these papers, T e (0) ~ 0.3 to 0.4 keV; these val ues are lower 
than the ~ 1 keV found by Cen fc Ostriker (1999) using 



hydrodynamical simulations with feedback effects. As a 
compromise between these results, we take T e (0) = 0.5keV 
and bg = 4, which corresponds to 



6^(0) = 0.0039. 



(19) 



Note that this is a fac tor of 2 lower t han used in Goldberg 
& Spjrgel] (1999) and |Cooray fc Hu| (1999). 

Needless to say, the resulting predictions should be 
taken as order-of-magnitude estimates only. As simula- 
tions improve, one can expect better values for the bias 
today and a more detailed modeling of its redshift and 
perhaps even scale dependence. 

3.2. Multipole Moments 

The multipole moments of the SZ effect under these sim- 
plifying assumptions can then be expressed a s a weighted 
projection of the density field ( Cooray fc Hu 1999): 



aZ(0) = J dhYr(n)e sz (h,0) 
d 3 k 



2tt 2 



£(k, n )/fWr*(k), 



(20) 



where 



7T 1 



I^(k)^W^(n)^--Fi(k), 

W sz (r) = -2^(r)f, 
in the Limber approximation and ( |Hrj 2000a) 

F l = {i + n K HlP/e)-^. 



(21) 



(22) 



The quantities take on a simple forms for a flat universe: 
ri — > l/k and Fi(k) — > 1. The Limber approximation 
breaks down for / <> 50 but is sufficient for our purposes. 

3.3. Power Spectrum 

The power spectrum of the SZ effect in this simplified 
model follows from equation ( pp| ) , 

cf z = ^/** > flt(*;n)[if z (*)] a , 



dr J Ps(l/d A ;r), 

d A 



(23) 



In the second line we have transformed the integration 
variable under the Limber correspondence: k — 1/cLa and 



dk 



Ft 



dr 
dA 



(24) 



We see that to go from the fiat to curved cosmologies in 
the Limber approximation one simply replaces the radial 
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Fig. 4. — top: SZ power spectrum from simulations compared to 
analytical predictions based on lin ear perturbation th eory (PT) and 
the non-linear scaling relations of Peacock fc Dodds| (1996; PD96). 
bottom: Errors on the binned power spectrum estimators for a single 
6° X 6° field; for a given experiment the errors should be scaled by 

~ 0.03/ sky . The sampling errors in the simulations is nearly equal 
to those of a Gaussian random field with the same power spectrum. 
The total noise including residual foregrounds and detector noise is 
also given for Planck. 



distance with the angular diameter distance in the inte- 
grand. 

In evaluating the SZ power spectrum, we have extended 
the SZ model to the non-linear regime by using the scal- 
ing for mulae for the nonl inear dark matter power spec- 
trum of Peacock fc Doddsj (1996). However, modeling the 



SZ effect with a scale-independent bias fa ctor will clearly 



break down deep in the non-linear regime. Refregier et al 



(1999) have shown that it is a reasonable approximation 
in the weakly non-linear regime (overdensities <^ 10) for 
z < 1 but can be in serious error outside of this range. 
As the weakly non-linear regime is the one of interest for 
anisotropies at I < 1000, we will use this approximation 
to test the effects of non-linearities. The predicted power 
spectrum in our fiducial model is shown in Fig. ^. 

3.4. Bispectrum 

The bispectrum of the SZ effect also follows from ex- 
pression ( |20|) 

r>mim 2 ni3 _ / \ 



(2 7 r) 3 <5 D (k 1 +k 2 +k 3 )B s (k u k 2 ,k 3 ) 



Here the density bispectrum should be understood as aris- 
ing from the full unequal time correlator 



(5(k 1 ;r 1 )^(k 2 ;r 2 )5(k 3 ;r3)) , 



(25) 



7 



where the temporal coordinate, which we temporarily sup- 
press, is evaluated in the Limber approximation (22). 

To further simplify this expression, we expand the delta 
function 



k 3 ) 



1 



,i(k 1 +k 2 +k 3 )-fir d 3 a . ^ ^ 2 gj 



(2tt) 3 

and employ the Rayleigh expansion 

gik-flr = 4vr ^^(fc r )^™*(k)^™(n) . (27) 

im 

We have assumed here a flat universe to simplify the 
derivation; as we have seen in the last section, we can 
promote the final result to a curved universe by replacing 
radial distances with angular diameter distances. 

With these relations, the angular integral over the di- 
rections of k 3 collapse to give 



By 



hi,. 



2 dr 



3=1 



xB(fc 1 ,fc 2 ,fc 3 )G™ m3 . 
where the Gaunt integral is 



G 



mim2m3 
I1I2I3 



dnY, mi Y, m2 Y l m3 

h *2 t3 



(28) 



(29) 



(2ii + l)(2i 2 + l)(2/ 3 + l) 

'1 h h \ I h h h 
/ 1 mi m 2 m.3 



Here, the quantities in parentheses are the Wigner-3j sym- 
bols whose p roperties are described in Appendix A of 
Cooray & Hu (1999). The integrals over the Bessel func- 
tions can again be done in the Limber approximation leav- 



ing 



mi 7712^3 



jjmim 2 m;i q 



[W**(r)]» A h h 



Note that only equal time contributions contribute in the 
Limber approximation. 

We can promote this result to a curved universe by re- 
placing radial distances with angular diameter distances 



d\ dA d A dA 



Finally, we can introduce the angular averaged bispectrum 
as 



E 



B hhh = 

mim2i)i3 

to obtain the final result 



h h h 
mi m 2 7713 



llhls 



(30) 



(2Zi + 1)(2Z 2 + 1)(2Z 3 + 1) (h l 2 h 





x dr 



4tt 

[W sz (r)f h l 2 I 



4 



*it't>i ir) - (31) 



One can alternately derive this relation by taking a flat- 
sky approach and using the general relation betw een the 
flat-sky and all-sky bispectra (see Appendix C, Hu 2000b). 

Equation (|l|) gives the SZ angular bispectrum in terms 
of the underlying density bispectrum. In second order per- 
turbation theory, the density bispectrum is in turn given 

by 



Bs{h, k 2 ,k 3 ;r) = F 2 (k 1: k 2 )P 5 (k 1 ;r)P s (k 2 - 1 r) 
+5 perm. , 



where 



i?n, v ^ 5 , k x • k 2 , 2 (k x • k 2 ) 

F 2 (ki,k 2 ) = - + — — + - p p 



(32) 



(33) 



Unfortunately, there exists no accurate fitting formula 
for the bispectrum of the density field in th e m ildly non- 
linear regime; we will employ simulations in §3.6 to address 
this regime. In the deeply non-linear regime, the density 
field obeys the hierarchical ansatz 



Bs(ki,k 2 ,k 3 ;r) 



Q 3 



[P(ki ;r)P(k 2 ;r)+5 perm.] , (34) 



where the power spectra are give n by the non-linear scal- 
ing of Peacock fc Dodds| (1996). Scoccimarro fc Frieman 
(1999) find that for power law power spectrum 



? 8 (n) = [4-2 n ]/[l + 2 



n+ll 



(35) 



Hui| (1999) suggests that for a general power spectrum one 



should replace n with the local linear power spectral index 
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Fig. 5. — top: Skewness in the simulations compared with second 
order perturbation theory (PT) and hyper-extended perturbation 
theory (HEPT). The smoothing is performed with an angular tophat 
of radius a. bottom: Errors on the skewness measurement for a 
single 6° X 6° field due to sampling errors and residual noise from 
Planck. 
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3.5. Skewness 

The simplest aspect of the bispectrum that can be mea- 
sured is the third moment of the map smoothed on some 
scale with a window W{a) 



I1I2I3 



(2Zi + l)(2Z 2 + l)(2I 3 + l) 



h h h 




4-7. 



B lll2l3 W h (a)Wi 2 (a)W h (4m 



where Wi are the multipole moments (or Fourier transform 
in a flat-sky approximation) of the window. For simplicity, 
we will choose windows which are either top hats in real 
or multipole space. 



It i* useful to define the skewiieaM parameter 



where the second moment is that of the SZ signal 



(6 2 (n;a)> 



47T t—' 

I 



(21 



l)CfWf(<7) 



(37) 



(38) 



The skewness in our fiducial model is shown for both the 
perturbation theory and HEPT predictions in Fig. |^. 

Since the density bispectrum in both the perturbative 
and non-linear regime scale as [^(/c)] 2 , the amplitude of 
the underlying density fluctuations roughly scale out of S3. 
However, the pressure bias does not: S3 cx 6" 1 . S3 thus 
provides an observable handle on the bias. This general 
point applies even if the bias is non-linear althoug h its 
interpretation will be no t be as straightforw ard (see Fry 



fc Gaztanaga 1993 and Mo, Jing fc White 1997 for its 



application in galaxy biasing) 

3.6. Numerical Simulations 

Since we are interested in the properties of the SZ effect 
in the weakly-nonlinear regime, cosmological simulations 
are required to recover the complete statistical properties 
of the signal and calibrate semi-analytic approaches for its 




- 

1 




Fig. 6. — One of 500 simulations of the SZ effect in the 
ACDM model for a 6° X 6° field-of-view. The range of the map is 
— lOOfiK , 25fiK with an rms of 9fiK and has an approximate angular 
resolution of 2'. Note the lack of obvious filamentary structures. 



low-order statistics. The simplified SZ model employed in 
this paper has the virtue that it is easy to simulate as it 
requires only dark matter and not the gas to model. Its 
main drawback of course is that results must be taken with 
a grain of salt due to missing physics. 

The realism of the basic approach can be improved by 
better calibrating the bias model against hydrodynamic 
simulations. One can envision going beyond the simple 
redshift dependent bias approach taken here to include 
scale dependence and stochasticity. Even accounting for 
these additional complications, simple dark matter simula- 
tions can continue to complement full hydrodynamic sim- 
ulations. Hydrodynamic simulations will always be more 
limited in dynamic range and sampling volume. Indeed, 
the current state of the art is limited a handful of realiza- 



tions across one or der of magn itude in physical scale (Re- 



frcgier et al. 1999; Beljak et al. 2000). A single simulation 



is then "stacked" on the linc-of-sight. Given the range of 
rcdshifts at which the SZ effect contributes, the simula- 
tion volume is traced many times for each line-of-sight. 
Moreover, the angular resolution decreases monotonically 
as one approaches the origin at z = 0. 

The reduction in dynamic range due to the angular pro- 
jection is a serious but not unfamiliar problem in cosmol- 
ogy. It occurs whenever t he kernel for t he projection spans 
cosmological distances. White fc Hu| (1999) introduced 
a technique of tiling multiple particle-mesh simulations 
which telescope along the line of sight to maintain a fixed 
angular resolution for the analogous problem in weak lens- 
ing. This also avoids the problem of over-r epresenting the 



filamentary structure of the map noted by Refregier et al 
(1999) 



We refer the reader to White fc Hu (1999) for details 



of the approach and tests of the method. The simulation 
all have a 256 3 mesh with 256 2 lines of sight for the ray 
tracing on a 6° x 6° field. Other relevant parameters are 
given in Tab. Ep the box size Lbox, the number of particles 
-^Vpartj the number of simulations run N s i m , the number of 
tiles of the given box size used iVtiic, the maximum redshift 
to which a given box is used, and the angular resolution of 
the mesh the maximum and minimum redshift used # mes h- 
Note that we cannot shrink the box size along the line-of- 
sight indefinitely since the fundamental mode of the box 
must be in the linear regime to provide accurate evolution. 
This implies that we lose angular resolution near the origin 
where a fixed physical scale subtends a large angle on the 
sky. Furthermore at the higher redshift the number of 
particles must be increased to eliminate shot noise from 
the initial conditions. Nonetheless, the tiling technique 
does a good job of maintaining angular resolution at all 
but the lowest redshifts. 

We construct 500 SZ maps from random combinations 
of the tiles in Tab. 1 for our statistical analysis; one real- 
ization is shown in Fig. |(| The average power spectrum is 
shown in Fig. ^ (top panel) and compared with the linear 
perturbati on theory predictio n and the non-linear scaling 
relation of Peacock & Dodds (1996). We have tested that 
the deficit of power at the low multipoles is an artifact of 
the finite field-of-view through monte-carlo realizations of 
the predicted power spectrum. The roll-off at high mul- 
tipoles is due to the spatial resolution in the simulations. 
This also explains the ~ 10% deficit at intermediate scales 
which comes from highly non-linear structure close to the 



Fig. 7. — Recovery of the SZ signal with Planck: left to right, modal SZ signal, signal + noise from primary anisotropics and foregrounds, 
and final recovered map from Planck. The signal map is that of Fig. H smoothed with a top- hot of radius 20'. 
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NOTES. — Numerical simulations in our ACDM cosmological model; 
see text for a description of these quantities. 



origin. Agreement is restored if one eliminates contribu- 
tions from overdensities > 10 in the predictions. Since our 
SZ model is at best valid in the weakly non-linear regime, 
these contributions should not be included anyway. 

Fig. H (top panel) shows the results for the skewness in 
the simulations compared with the second order perturba- 
tion theory and HEPT predictions. The agreement here 
is worse, but is still sufficient for our purposes, given the 
crudeness of the underlying model for the SZ effect itself. 

We can address sample variance questions from the scat- 
ter of the results in the individual realizations. Sampling 
errors for the power spectrum and skewness are shown in 
the bottom panels of Fig. ^ and [| respectively. Since these 
are for individual 6° x 6° planes, they should be scaled by 

— 1/2 

~ 0.03/ sky for a given experiment. Sampling errors are 
one source of noise that we will include in the signal-to- 
noise calculations in the next section. 

4. ESTIMATING THE SIGNAL-TO-NOISE 

With the SZ signal estimated from the simple bias model 
of §|| and residual noise calculated from the foreground 
model and subtraction techniques of §|^, we can now es- 
timate the signal-to-noise for the detection of the SZ ef- 
fect. In Fig. 0, we illustrate the foreground subtraction 
technique on simulated Planck maps. The signal-to-noise 
in the maps is of order one for features spanning tens of 
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Fig. 8. — Cumulative signal-to-noise in the measurement of the 
SZ power spectrum with Boomerang, MAP and Planck as a func- 
tion of maximum I. The solid line is the maximum signal-to-noise 
achievable in a perfect experiment (see text). 



arcminutes. We shall here show that this level of signal-to- 
noise is more than sufficient for the purpose of extracting 
measurements of the low order statistics of the SZ signal. 

4.1. Power Spectrum 

The signal-to-noise in the power spectrum per multipole 
(l,m) mode is simply 



1 /cp 

2 \Cf ot 



(39) 



Here, Cf ot is the power spectrum of all contributions in 
the SZ map, 



C, sz 



N, 



(40) 



where recall that the residual noise Ni was defined in equa- 
tion (|l0|) and includes contributions from detector noise. 

Assuming Gaussian statistics for the signal and noise, 
each mode is independent so that the total signal-to-noise 



10 



is the quadrature sum 



1000 



-) =^*>> ' 



fCf^ 



(41) 



This quantity gives the variance of the total power mea- 
surement in the SZ effect, including sample variance. 
S/N 3> 1 means that one has a precise measurement of the 
power spectrum not simply a highly significant detection. 
In Fig. ||, for the Boomerang, MAP and Planck experi- 
ments as a function of the maximum I mode included in 
the sum. We also show the ultimate limit of perfect fore- 
ground and noise removal where Cf ot = Cf z and / s k y = 1. 
We will refer to this case here and below as a "perfect 
experiment" . 

With our fiducial choice of the gas bias, Planck should 
have a highly significant detection of the total signal. One 
should bear in mind that the bias parameter is still 
highly uncertain and that the S/N scales as b\. Never- 
theless even a relatively large reduction in the average gas 
temperature or density bias will not make the signal unde- 
tectable in principle. In practice, however remember that 
one is then relying on a precise subtraction of the noise 
bias in the measurement of C* ot , which in turn requires 
that the power spectrum of the dust and other residual 
foregrounds lurking at least at the 10% level in rms (1% 
in power) are determined comparably precisely. 

If the fiducial SZ bias is close to correct, the high total 
single-to-noise in Planck can be used to break the mea- 
surement into bands in I and recover the band power with 
errors 
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(42) 
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We give an example from monte carlo realizations of the 
Gaussian noise and sample variance from the simulations 
in Fig. H Note that these are errors for a 6° x 6° section 
of the sky and should be scaled by 
Planck. 

These signal-to-noise estimates assume that both the 
signal and noise are Gaussian. Of course in reality the 
SZ signal is non-Gaussian. In general, gravitational col- 
lapse correlates the amount of power in density fluctu- 
ations across all scales in the non-linear regime. How- 
ever since the SZ effect probes many independent density 
fluctuations along the line-of-sight, the central limit theo- 
rem ensures that the SZ signal is far more Gaussian than 
the density field. We can test how much this affects the 
signal-to-noise with our simulations. Shown in Fig. |J is 
the sampling errors on the band powers from the simula- 
tions themselves as compared with those from Gaussian 
realizations of the same power spectrum. The excess vari- 
ance over the Gaussian limit is small on the relevant scales 
given detector noise limitations from Planck. 

4.2. Skewness 

The overall signal-to-noise for the measurement of the 
third moment of SZ effect is 
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(43) 
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Fig. 9. — Cumulative signal-to-noise in the measurement of the 
third moment, 3 , with top hat smoothing in multipole space (i.e. 
truncation above i ma x)- The HEPT approximation to the bispec- 
trum is assumed here. MAP and Boomerang (not shown) have 
signal-to-noise values less than 0.1 everywhere. 



where the variance is given by 



Var = 



1 



(2Ji + l)(2l a + l)(2l3 + l) 
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h h h 
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In Fig. 0, we show the signal-to-noise for a measurement 
of the third moment as calculated under the HEPT. We 
compare the signal-to-noise in Planck with the ideal case 
of perfect removal of foregrounds and detector noise, and 
full sky coverage. We use here a tophat window in mul- 
tipole space out to ^ max to conform with other signal-to- 
noise considerations. Cosmic variance and Planck detec- 
tor noise reduces the signal-to-noise values both at the low 
and high end for i max values respectively. For Planck, the I 
values in the range of few hundred to ~ 1000 provides the 
maximal signal-to-noise for a measurement of the skew- 
ness. This corresponds to smoothing scales a ~ 10'-30' for 
tophat windows in angular space (c.f. Fig. g). For MAP 
and Boomerang, the signal-to- noise values are <; 0.1, sug- 
gesting that a detection of SZ skewness is not likely to be 
possible in these two experiments. 

Again equation (|44|) assumes Gaussian statistics for the 
variance and ignores the sample variance of the third mo- 
ment itself. We test this approximation in Fig. H and find 
that it is reasonable given the level of residual noise for 
Planck. In constructing an estimator for S3, it is im- 
portant to remove the noise bias since noise variance will 
always reduce the skewness in the map. We do this by 
multiplying the estimator by ((©t ot ) / (©sz)) 2 - 

Finally, note that in the noise-dominated regime the 
signal-to-noise in S3 scales strongly with the gas bias 
S/N oc b\ , so that the detectability of this effect depends 
strongly on currently uncertain assumptions. 

4.3. Bispectrum 

The full bispectrum of the SZ effect contains all of the 
information about its three-point correlations induced by 
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Sperge] 1999; Cooray fc Hu 1999). Again the cosmic vari- 
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Fig. 10. — Cumulative signal-to-noise for the detection of SZ bis- 
pectrum as a function of £3 multipole. The solid line is the maximum 
signal- to- noise achievable in a perfect experiment. 



the growth of structure beyond the linear approximation. 
The skewness is simply one, easily measured, aspect of the 
bispectrum. The full signal-to-noise ratio of the bispec- 
trum is 



/sky ^2 
li ,l2,h 



ll 12 13 



(45) 



where C* ot follows Eq. (40). We plot the bispectrum cu- 
mulative signal-to- noise as a function of signal ^3, su mmed 
over li and li- We refer the reader to Cooray fc Hu| (1999) 



for a detailed discussion on the bispectrum, its variance 
and the calculation of signal-to-noise ratio. 

In Fig. |l^, we show the expected cumulative signal- 
to-noise for the SZ bispectrum in Boomerang, MAP and 
Planck data and a perfect experiment. The signal-to- noise 
is calculated under the HEPT approximation for the un- 
derlying density field. As shown, MAP and Boomerang 
allow reasonable limits to be placed on any non-Gaussian 
signal in the SZ effect while Planck allows a strong possi- 
bility for a detection. 

Again the same caveats as to the sensitivity of the S/N 
estimate to the underlying assumptions that applied for 
the skewness also apply here. Moreover, measuring all the 



config urations of the bispectrum will be a formidable com- 
putational challenge as will control over systematic etlects 
in the experiments. 

4.4. Lensing Correlation 

The SZ effect and weak gravitational lensing of the CMB 
both trace large-scale structure in the underlying density 
field. By measuring the correlation, one can directly test 
the manner in which gas pressure fluctuations trace the 
dark matter density fluctuations. The correlation vanishes 
in the two-point functions since the lensing does not affect 
an isotropic CMB due to conservation of surface bright- 
ness. 

The correlation manifests itself as a non-v anishing bis- 
pectrum in the CMB at RJ frequencies (Goldberg & 



ance from the primary anisotropics presents an obstacle 
for detection of the effect above the several arcminute scale 
(/ ~ 2000). With the multifrequency cleaning of the SZ 
map presented here one can enhance the detectability of 
the effect. 

Consider the bispectrum composed of one ctim from the 
cleaned SZ map and the other two from the CMB maps. 
Call this the SZ-CMB-CMB bispectrum. The noise vari- 
ance of this term will be reduced by a factor of Cf ot /C ; CMB 
compared with the CMB-CMB-CMB bispectrum. As one 
can see from Fig. || this can be up to a factor of 10 3 in the 
variance. Details for the calc ulation of the CMB-CMB- 
CMB bispectrum are given in |Cooray fc Hvj (1999). Here 
we have updated the normalization for SZ effect, taken 
/sky = 0.65 for Planck's useful sky coverage, and compared 
the S/N of the two bispectra. As shown, the measurement 
using foreground cleaned Planck SZ and CMB maps has 
a substantially higher signal-to-noise than that from using 
the Planck CMB map alone for multipoles I ~ 10 2 — 10 3 . 

Beyond the improvement in signal-to-noise, however, 
there is an important advantage in constructing the SZ- 
lensing bispectrum using SZ and CMB maps. A mere 
measurement of the bispectrum in CMB data can lead 
to simultaneous detection of non-Gaussianities through 
processes ot her than just SZ-lcns ing cross-correlation. As 
discussed in Gloldberg fc Spergej (1999) and extended in 



Cooray & Hu (1999), gravitational lensing also correlates 



with other late time secondary ani sotropy contrib utors 
such as integrated Sachs- Wolfe (ISW; |Sachs fe Wolfe] 1967) 
effect and the reionized Doppler effect. In addition to lens- 
ing correlations, non-Gaussianities can also be generated 
throu gh rcionization and non-li near growth of pert urba- 
tions (fepcrgel fc Goldberg 1999; |Goldberg fc Spergel| 1999; 



Cooray & Hu 1999). Bispectrum measurements at a signle 



frequency can result in a confusion as to th e relative con- 
tribution from each of these scenarios 



In 



Cooray & Hr 



(1999), we suggested the possibility of using differences in 
individual bispectra as a function of multipoles, however, 
such a separation can be problem atic given that these dif- 
ferences are subtle (e.g., Fig 6 of Cooray & Hu 1999). 

The construction of the SZ-lensing bispectrum using SZ 
and CMB maps has the advantage that one eliminates all 
possibilities, other than SZ, that result in a bispectrum. 
For effects related to SZ, the cross-correlation of lensing 
an d SZ should p roduce the dominant signal; as shown 
in Cooray fc Hu| (1999), bispectra signal through SZ and 



rcionization effect s, such as Ostriker-Vishniac (OV; Os 



triker & Vishniac 1986), are considerably smaller. 



Conversely, multifrequency cleaning also eliminates the 
SZ contribution from the CMB maps and hence a main 
contaminant of the CMB-CMB-CMB bispectrum. This 
assists in the detection of smaller signals such as the ISW- 
lensing correlation, Doppler-lensing correlation or the non- 
Gaussianity of the initial conditions. Such an approach is 
highly desirable and Planck will allow such detailed studies 
to be carried out. 

A potential caveat is that as noted above, the full bis- 
pectrum in an all-sky satellite exp eriment will be difficult 
to measure. Zaldarriaga fc Seljak| (1999) have developed a 



reduced set of three-point statistics optimized for lensing 
studies, based on a two point reconstruction of the lensing- 
convergence maps from temperature gradient information. 
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They show that most of the information is retained in these 
statistics. Multifrequency cleaning improves the signal-to- 
noise for these statistics by exactly the same factor as for 
the full bispectrum. 

5. DISCUSSION 

We have studied the prospects for extracting the statis- 
tical properties of the Sunyaev-Zel'dovich (SZ) effect asso- 
ciated with hot gas in large-scale structure using upcoming 
multifrequency CMB experiments. This gas currently re- 
mains undetected but may comprise a substantial fraction 
of the present day baryons. The SZ effect has a distinct 
spectral dependence with a null at a frequency of ~ 217 
GHz compared with true temperature anisotropies. This 
frequency dependence is what allows for effective sepa- 
ration of the SZ contribution with multifrequency CMB 
measurements . 

As examples, we have employed the frequency and noise 
specifications of the Boomerang, MAP, Planck experi- 
ments. The MAP satellite only covers frequencies at 
RJ part of the frequency spectrum. Consequently, only 
Boomerang and Planck can take full advantage of multifre- 
quency separation of the SZ and primary anisotropies. We 
have evaluated the detection threshold for SZ power spec- 
trum measurements (see Fig. |). Boomerang and MAP 
should provide limits on the degree scale fluctuations at 
the several /xK level in rms; Planck should be able to de- 
tect sub fiK signals. 

The expected level of the SZ signal in the fiducial ACDM 
model is still somewhat uncertain. We have employed a 
simple bias model for the pressure fluctuations , roughly 



norma lized t o recent hydro dynamic simulations (Refregiei 
et al. 1999; Seljak et al. 2000), and calculated the re 



suiting signal using analytic scaling relations and particle- 
mesh dark matter simulations. As hydrodynamic simula- 
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Fig. 11. — Cumulative signal-to- noise in the measurement of 
the SZ-weak gravitational lensing cross-correlation through the bis- 
pectrum measurement in CMB data. Compared are the expected 
signal-to-noise with (SZ-CMB-CMB) and without (CMB-CMB- 
CMB) multifrequency isolation of the SZ effect for Planck and a per- 
fect/cosmic variance limited experiment. Multifrequency isolation 
provides additional signal-to-noise and the opportunity to uniquely 
identify the bispectrum contribution with the SZ effect. 



tions improve, these techniques can be extended with more 
sophisticated modeling of the bias. They complement hy- 
drodynamic simulations by extending the dynamic range 
and simulated volume, the latter being important for ques- 
tions of sample variance. 

Assuming this simplified model of the SZ signal, Planck 
should have signal-to-noise per multipole of order unity for 
I < 1000. Although the recovered maps are then somewhat 
noisy, they are sufficient for precise determinations of low 
order statistics such as the SZ power spectrum, bispectrum 
and skewness (see Figs. f|-[l(]). The skewness in principle 
can be used to separate the pressure bias from the under- 
lying amplitude of the density fluctuations. The full bis- 
pectrum contains significantly more information but will 
be difficult to extract in its entirety. Current methods for 
measuring the bispectrum, tested with the COBE data, 
have concentrated a t measuring sp ecific modes such as 
l\ = l 2 = h = I (Ferreira et al. 1998). More work 



will clearly be required, especially in understanding the 
systematic errors at a sufficient level, but the wealth of 
information potentially present in the bispectrum should 
motivate efforts. 

Note however that the non-Gaussianity in the SZ signal 
is not very strong due to the fact that it is constructed 
from many independent pressure fluctuations along the 
line of sight. As a consequence, we expect that signal-to- 
noise ratios can be estimated by Gaussian approximations, 
but that techniques that try to impr ove the SZ-pr imary 
separation based o n non-Gaussianity ( Hobson et al| 1998; 
Aghanim fc Forni| 1999) may not be particularly effective 



for this signal. 

We caution the reader that our oversimplification of the 
SZ signal can cause problems for a naive inter pretation of 
future detections. For example, |Seljak et al. (2000) find 
that the SZ power spectrum in their simulations is domi- 
nated by shot noise from the rare hot clusters not included 
in our modeling. Fortunately since these contributions are 
highly non-Gaussian, they can can readily be identified 
and removed. At the very least, X-ray bright clusters 
can be externally identified and removed; this has been 
shown to substantially reduce the shot noise contribution 
(Komatsu & Kitayama 1999). The effect we are modeling 
should be understood as the signal in fields without such 
clusters. 

Another means of separating the SZ signal from large- 
scale structure from that of massive clusters is to cross 
correlate it with other tracers of large-scale structure that 
are less sensitive to highly overdense regions. An added 
benefit is that such a cross-correlation will also empiri- 
cally measure the extent to which pressure fluctuations fol- 
low mass fluctuations. The CMB anisotropies themselves 
carry one such tracer in the form of the convergence from 
weak lensing. It ma nifests itself as a thr ee-point correla- 
tion or bispectrum (Goldberg & Spergel 1999) but with- 
out frequency information it is severely sample-variance 
limited due to confusion noise from primary anisotropies. 
Measuring the SZ-lensing correlation using the cleaned SZ 
maps improves the signal-to-noise for the detection by over 
an order of magnitude at degree scales. Furth ermore, the 
techniques introduced by Zaldarriaga & Seljak (1999) pro- 



vide a concrete algorithm for extracting most of the three- 
point signal without recourse to measuring all the config- 
urations of the bispectrum. Conversely, SZ removal from 
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the CMB maps themselves can assist in the detection of 
other smaller bispectrum signals by eliminating one source 
of confusion noise. 

The cross-correlation coefficient between the SZ effect 
and CMB w eak lensing is relatively modest (~ 0.5, see 



Seljak et al| 2000). This is due to the fact that the SZ 
effect is a tracer of the nearby universe while CMB lens- 
ing is maximally sensitive to structure at z ~ 3. A higher 
correlation is expected if SZ is cross-co rrelated with an 
external probe of low redshift structure. Peiris fc Spergej 
(2000) suggested the cross-correlation of MAP CMB data 
and Sloang galaxy data. An improved approach would be 
to use the Planck derived SZ map rather than a CMB 
map. Using a SZ map reduces noise from the primary 
anisotropies and guarantees that any detection is due to 
co rrelations with th e SZ effect. Extending the calculations 



Peiris & Spergel (2000) with the Planck generated SZ 



map, we find signal-to- noise ratios which are on average 
greater by a factor of ~ 10 when compared to signal-to- 
noise values using MAP CMB map. In fact with redshifts 
for galaxies, Planck SZ map can be cross-correlated in red- 
shifts bins to study the redshift evolution of the gas. Other 
promising possibilities include cross correlation with soft 
X-ray background measurements, as well as ultraviolet and 
soft X-ray absorption line studies. 

All these considerations imply a bright future for SZ 
studies of the hot gas associated with large-scale structure 
with wide-field multifrequency CMB observations. Its de- 
tailed properties should be revealed in its non-Gaussianity 
and correlation with other tracers of large-scale structure. 
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